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Abstract 

In systems where deposition rates are high compared to diffusion, desorption and other mecha- 
nisms that generate correlations, a crossover from random to correlated growth of surface roughness 
is expected at a characteristic time to- This crossover is analyzed in lattice models via scaling ar- 
guments, with support from simulation results presented here and in other authors works. We 
argue that the amplitudes of the saturation roughness and of the saturation time tx scale as 
tQ-/"^ and respectively. For models with lateral aggregation, which typically are in the Kardar- 
Parisi-Zhang (KPZ) class, we show that to ~ P^^i where p is the probability of the correlated 
aggregation mechanism to take place. However, to ~ p "^ is obtained in solid-on-solid models with 
single particle deposition attempts. This group includes models in various universality classes, 
with numerical examples being provided in the Edwards- Wilkinson (EW), KPZ and Villain-Lai- 
Das Sarma (nonlinear molecular-beam epitaxy) classes. Most applications are for two-component 
models in which random deposition, with probability 1 — p, competes with a correlated aggregation 
process with probability p. However, our approach can be extended to other systems with the same 
crossover, such as the generalized restricted solid-on-solid model with maximum height difference 
S, for large S. Moreover, the scaling approach applies to all dimensions. In the particular case of 
one-dimensional KPZ processes with this crossover, we show that to ~ i'^^ and u ^ A^/^, where v 
and A are the coefficients of the linear and nonlinear terms of the associated KPZ equations. The 
applicability of previous results on models in the EW and KPZ classes is discussed. 

PACS numbers: PACS numbers: 81.15Aa, 68.55. Jk, 05.40.-a, 05.50.+q 



I. INTRODUCTION 



The large number of technological applications of thin films and multilayers motivates the 
study of continuous and atomistic models for the growth of those structures . Their 

morphology is often the product of a competition between different growth dynamics, thus 
theoretical models representing such features received increasing attention in recent years. 
Some examples are the models with aggregation of different species of particles la l a or 
those mixing different microscopic aggregation rules for the same species 
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li. 



In processes starting from a fiat surface where the deposition rate is high compared 

to ad-particle diffusion coefficients, a random growth is initially observed, with negligible 

spatial correlations in the local thicknesses. Subsequently, diffusion, desorption and other 

mechanisms introduce surface correlations and, consequently, a crossover from random to 

correlated growth is observed. The simplest quantitative characteristic of the film surface 

which reveals this crossover is the average roughness (or interface width), defined as the 

r/. -\2\1^/^ 

rms fluctuation of the height h around its average position h: W{L,t) = (Ih — hj) 
The aim of this work is to study scaling relations for the surface roughness in lattice models 
which present this crossover. 

In a random, uncorrelated growth, the roughness increases as 

Wr ~ t^/'. (1) 

On the other hand, in correlated growth processes, it is expected to obey the Family- Vicsek 
scaling relation [l^ 

W{L,t) ^ Aiy (j-y (2) 

where L is the system size, a is the roughness exponent, A is a model-dependent constant, 
/ is a scaling function such that / ~ 1 in the regime of roughness saturation [t oo) and 
tx is the characteristic time of crossover to saturation, tx scales with the system size as 

tx ^ BU, (3) 

where z is the dynamic exponent and B is another model-dependent amplitude. For t <^tx 
(after a possible crossover to correlated growth), the roughness scales as 

W ^Ct^, (4) 



where C is constant and j3 = a/z is the growth exponent. In this growth regime, /(x) ~ 
in Eq. (0). 

Fig. 1 illustrates the time evolution of the surface roughness in systems which present 
a crossover from random to correlated growth at a characteristic time tn (this time is also 
called ... a.d b. otKe. autKo. QSQ). Fo. . < W ...ease. a. E,, ffl. and the 
reduced slope in the subsequent regime (to ^ ^ ^ ^x) is a signature of the smoothing effect 
of correlations. In systems where a correlation mechanism and a random growth mechanism 
are simultaneously present, a different balance of these mechanisms changes the crossover 
time to, which affects the amplitudes A, B and C (Eqs. El El and Ej) . 

This was observed by Albano and co-workers p, [lo|, lllj , who studied numerically two- 
component models in which different rules for the aggregation of the same species are chosen. 



with coinplementarv probabilities (they are simply called competitive models by several 
authors fl H H Q Q, Q, Q). The first component of those models, with probability 
1 — p, was random deposition (RD): particles attach to the top of the column of incidence 
independently of the neighboring heights. In their first competitive model (called RD-BD 
model), the second component was ballistic deposition (BD) with probability p. In 
their second competitive model (called RD-RDSR), the correlated component was random 
deposition with surface relaxation (RDSR) [21]. As expected from the absence of correlations 
in RD, the universal scaling exponents a, (3 and z of the competitive models are those of 
the universality class of the correlated component However, the scaling amplitudes are 
affected, and scale with p as 

A-p-^ (5) 

B ~ p-y, (6) 

and 

C ^p-^. (7) 

or the RD-BD model, 5 = 1/2 and y = 1 were obtained in all substrate dimensions d 
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11|. For the RD-RDSR model, 5 = 1 and ?/ = 2 were obtained Q- In both 
competitive models, the exponent 7 depended on d. 

In a recent work, Braunstein and Lam explained the differences between those systems 
through scaling arguments which account for the average height increase and fluctuations 
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during the time interval between two correlated deposition events. Previou sly, the equation 



representing the RD-RDSR model in the continuum limit was also derived 



14| . Since those 



22 1 and in the Kardar- 



works analyzed particular models in the Edwards- Wilkinson (EW) 
Parisi-Zhang (KPZ) 



23| classes, the picture that emerged from those results was that the 
exponents 6 and y are related to the universality class of the dominant process: 6=1/2 
and y = 1 in the crossover to KPZ, 6 = 1 and y = 2 in the crossover to EW. The dominant 
behavior of another competitive model in Ref. |l6|] was not trivial to infer, but the crossover 
exponents to KPZ and EW were the same. However, recent numerical work by Kolakowska 
et al show that 6=1 and y = 2 is also found in two-component models in the KPZ 

class. Thus, a complete theoretical explanation of the values of these crossover exponents is 
still lacking. 

Here we will present a scaling approach which provides such explanation through a con- 
nection of roughness amplitudes, crossover times and the microscopic rules of the lattice 
models. It can be applied to all spatial dimensions and to all types of correlated growth, 
and agrees with our numerical results for models in three different universality classes and 
with numerical results from other authors works [lol llll. Ifoj]. 

We will argue that A ~ to^^^ and B in the models with random to correlated 

crossover, so that to acts as a time dilatation factor. The relations between to and the 
parameters of the discrete models are subsequently obtained by scaling arguments, partially 
rephrasing those of Braunstein and Lam . These arguments lead to a separation of the 
lattice models in two groups: the first one includes solid-on-solid models with single particle 
deposition attempts, for which 6 = 1 and y = 2, and the second one includes models with 
lateral aggregation, for which 5 = 1/2 and y = 1. This classification is independent of the 
universality class of the correlated process, consequently KPZ processes may be found in 
both groups, i. e. with different crossover exponents. Moreover, our theoretical approach 
comprises not only two-component models involving RD, but can be extended to other 
models with the same crossover. One example is the generalized restricted solid-on-solid 

n 

(RSOS) model, with maximum height difference S between the neighbors jl9(], which presents 
that crossover for large S {S^^ substitutes p in Eqs. El and [7j) and which will be studied 
numerically here. 

The theoretical analysis is motivated in Sec. II, with the presentation of several discrete 
models showing the random to correlated crossover and the discussion of their simulation 
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results. The scaling approach is presented in Sec. Ill, and in Sec. IV we summarize our 
conclusions. 



II. LATTICE MODELS AND SIMULATION RESULTS 

First we recall the models with crossover from random to correlated growth previously 
studied by other authors. 

The one . RD-BD Q. pure BD (p ^ 1), particle s«c.. at contact 
with a nearest neighbor, as illustrated in Fig. 2a, which leads to the formation of a porous 
deposit. Simulations of the RD-BD model gave 5^1/2 and yK.l'm.d = l,d = 2 and = 3 
10, Q|, while 7 depended on d. From the Family- Vicsek relation Q and the p-scaling for 



the amplitudes A, B and C (Eqs. El and [7j) , Albano and co-workers proposed that 

/5y-5 + 7 = (8) 



Here /? is the growth exponent 



in any rf, which agrees with their simulation results ^ 
of the KPZ class. 

The KPZ equation, which describes BD in the continuum limit, is 

^-^ = yV'h+^{Vhf + r^{x,t), (9) 

where h is the height at the position a; in a (i-dimensional substrate at time t, u represents a 
surface tension, A represents the excess velocity and r/ is a Gaussian noise P, 13] with zero 
mean and co-variance {i] {x,t) t]{x' ,t')) = D6'^{x — x')6 {t ~ t'). In d = 1, the exact KPZ 
exponents are a = 1/2, (3 = 1/3 and z = 3/2, and in d > 2 approximate values are given in 



Refs. 



2J and 



2^. In models with a crossover from random to KPZ growth, small values 
of u and A are expected in the corresponding KPZ equation. 

The second model studied by those authors was RD-RDSR [3]. In pure RDSR {p = 1), 
he incident particle diffuses to the column with minimum height in its nearest neighborhood 



21[. RDSR is described by the EW equation, which corresponds to Eq. El with A = 0. In 

were obtained for the RD-RDSR 

model, while the exponent 7 also depended on d. 

Other models involving competition with RD in c? = 1 were recently proposed in Refs. 



Q and |l7| 



One example is a model whose correlated component allows aggregation of 
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17l |. It is in the KPZ class, but has 5^1 and 



the incident particle only at surface minima 
1/ ~ 2. 

In the following, we will present our numerical results for three models devised to broaden 
the investigation on the crossover exponents. 

The first one is also a two-component model, involving RD with probability 1 — p. With 
probability p, the aggregation is possible only if the height of the column of incidence does 
not exceed the heights of any of its neighbors, otherwise the aggregation attempt is rejected, 
as shown in Fig. 2b. In other words, aggregation is possible only in valleys or plateaus. 
This model mimics a competition between RD and RSOS deposition, thus we call it the 
RD-RSOS model. We refer to the correlated component as RSOS because it works against 
the formation of large local slopes, as illustrated in Fig. 2b (there, the attempts at columns 
3 and 9 are rejected because one or more neighbors have smaller heights). In the pure RSOS 
model j3| {p = 1), the above aggregation rule implies Ahmax = 1 between neighboring 
columns. 

The RD-RSOS model was simulated in one- dimensional lattices with 32 < L < 512 for 
some values of p in the range 0.12 < p < 0.4. lO'' realizations were obtained for the smallest 
lattices and 10'^ for the largest ones. For the same values of p, the model was simulated in 
substrates with L = 8192 during the random and the KPZ growth regimes. 

In Fig. 3 we show the time evolution of the roughness for p = 0.15 and p = 0.3 in lattices 
with L = 256. For both values of p, it is clear that the roughness behaves as in the sketch of 
Fig. 1. The saturation regimes are clearer in the linear plot of the inset, which shows that 
the saturation roughness approximately doubles when p is reduced from 0.3 to 0.15. Using 
Eqs. El and this suggests 6^1. The parallel lines of the KPZ growth regimes illustrate 
the behavior described by Eqs. (jH) and (|7j), with the amplitude C increasing as p decreases. 

In order to obtain reliable estimates of the crossover exponents, we analyzed the effects 
of finite L, t and p, taking the limits L oo, t oo, p ^ when appropriate. This 
procedure was proved to be essential to avoid erroneous conclusions on the class of several 
growth models (see, e. g. the analysis of BD in j^])- For a fixed value of p, our first step 
is to extrapolate Wsat/L°' to L ^ oo, using a = 1/2, as illustrated in Fig. 4a for p = 0.25. 
From Eq. Q, the asymptotic value of that ratio is the amplitude A{p) {A ^ 3.18 in Fig. 
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4a). For each pair of subsequent probabilities and p", we define effective exponents 

_ In [A{p')IA{p")] 

Inp'Vp' ^ ^ 

where p is an average probabihty: 

V = \/(pyo- (11) 

As p' and p" decrease, p — and Sp 6. In Fig. 4b we show Sp versus p^, which suggests 

The first step to estimate 7 is to extrapolate W/ 1^^^ in the growth regime, using the data 
from lattices with L = 8192. For fixed p, that ratio converges to C{p) as t ^ 00 (Eq. |^. 
This is illustrated in Fig. 4c for p = 0.2, where we obtain C ~ 1.22 (this extrapolation 
represents the long time behavior in an infinitely large substrate). Effective exponents jp 
were calculated from C{p) along the same lines of the calculation of Sp in Eq. (fTTl|) . They 
are shown in Fig. 4d as a function of p^, which suggests 7~l/3asp^0. The estimates 
of 6 and 7 and Eq. (jH)) give y ^ 2 for the RD-RSOS model. 

Thus, although the RD-RSOS model is in the KPZ class, similarly to the RD-BD model, 
its crossover exponents 6 and y are the same of the RD-RDSR model, which is in the EW 
class. 

The second model analyzed here does not involve competition of ag greg ation rules. It 
is called generalized RSOS model, and was originally proposed in Refs. j2a] and 28 1. The 
incident particle can aggregate at a certain column only if the height differences between 
neighboring columns do not exceed an integer value S, otherwise the deposition attempt is 
rejected. The version with 5 = 3 is illustrated in Fig. 2c. For large S, with an initially 
flat substrate, random growth occurs until a significant fraction of neighboring columns has 
height difference S. Subsequently, KPZ growth is observed due to the rejection of deposition 
attempts. 

The generalized RSOS model was studied numerically by Chien et al who obtained 
to ~ S"^'^^ for large S, in agreement with their scaling arguments, which give to ~ S"^- 
However, the roughness amplitudes were not calculated there. 

Here it was simulated until saturation in one- dimensional lattices with 32 < L < 512, for 
several values of 5* in the range 4 < 5* < 32. The number of realizations was 10^ for the 
smallest lattices and 10^ for the largest ones. We also simulated the model in L = 8192 up 
to t ^ tj., with 10^ realizations for each S. 



For any 5", the height difference A/i = /i^+i — hi of neighboring columns can assume 
25 + 1 different values. However, deposition attempts are rejected only when Ah is —S or 
+S, because those attempts would lead to height differences — (S* + 1) and S" + 1. For large 
S, it is reasonable to assume that all values of Ah have nearly the same probability, thus the 
probability of rejecting aggregation is of order ^ S^^. Since aggregation rejection is 
the mechanism to spread correlations, plays the same role of p in the other competitive 
models. Thus we assume that 

A{S) ~ S\ (12) 
B{S) ~ S\ (13) 

and 

CiS) ~ S"' (14) 

in the generahzed RSOS model. 

Estimates of the crossover exponents 5 and 7 were obtained along the same lines of the 
RD-RSOS model described above. First, for fixed S, the extrapolation of Wsat/ L^^'^ provided 
estimates of the amplitude ^(5*). Finite-size estimates of the exponent 6 are given by 

ln[A{S)/A{S/2)] 

= • 

Their values are shown in Fig. 5a as a function of 1/S, suggesting 5^1 asymptotically. 
Amplitudes C{S) were obtained from the extrapolation of W/t^^^ in the KPZ growth regime 
of large substrates, and effective exponents 75- were defined analogously to Eq. (fT^ . They 
are shown in Fig. 5b as a function of 1/5*, suggesting 7^1/3 asymptotically. Those values 
also lead to ?/ ~ 2. 

The third model analyzed here is a two-component one which belongs to the class of 
the Villain-Lai-Das Sarma (VLDS) equation 29|, |30|. Again, RD has probability 1 — p. 
With probabihty p, aggregation is allowed only at valleys or plateaus, similarly to the RD- 
RSOS case. However, if aggregation is not possible at the column of incidence, then the 
incident particle migrates to the nearest column in which that condition is satisfied, and is 
irreversibly attached there. When p = 1, we obtain the conserved restricted solid-on-solid 
(CRSOS) model of Kim et al jsi)], where heights differences between neighboring columns 
do not exceed 1. However, with p < 1, differences of column heights larger than one appear 
due to the RD component. The competitive model will be called RD-CRSOS. 
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The original CRSOS model, as well as the RD-CRSOS, are represented in the continuum 
limit by the VLDS equation (also called nonlinear molecular-beam epitaxy equation) 3^ [s^ 

— = jj^V^h + KV\Vhf + r7(x, t), (16) 

where and A4 are constants. The best known estimates of scaling exponents for the VLDS 
class were obtained from extensions of the CRSOS model; in c? = 1, they are a = 0.93 and 
z = 2.88 I24I, which gives /? = 0.323. 

The RD-CRSOS model was simulated in one-dimensional lattices with the same values 
of p and L of the RD-RSOS model. In order to estimate the amplitudes A and C (Eqs. |21 
and|^, we used the above estimates of a and (3. In Fig. 5c we show 6p versus p^, which 
suggests 5 ~ 1, although the oscillations in the effective exponents do not allow a reliable 
extrapolation to p — > 0. In Fig. 5d we show 7^ versus p^, which suggests 7 ~ 1/3. These 
values also lead to y ~ 2. 

Thus, we have shown that a model in the VLDS class also has the crossover exponents of 
the previous solid-on-solid models in the EW and the KPZ classes, differing only from the 
value of the RD-BD model. 



III. SCALING THEORY FOR THE CROSSOVER FROM RANDOM TO CORRE- 
LATED GROWTH 

First we discuss the universal features of this crossover in lattice models in any substrate 
dimension d, despite the illustrations are all given in d = 1 for simplicity (e. g. in Fig. 2). 
We consider here models without additional crossovers between different growth dynamics. 

It is reasonable to assume that is the time in which the roughness of random growth 
(Eq. P) matches the growing roughness of the correlated process, Eq. 1^ (see Fig. 1). This 
gives 

C ~ to'/'-^. (17) 

Now to is the characteristic time for the onset of correlations among neighboring columns, 
which otherwise randomly grow. On the other hand, in a pure correlated model, the time 
6t = 1 for deposition of one monolayer is enough to produce such correlations. Consequently, 
in models with the random to correlated crossover, we expect that all characteristic times are 
scaled by a factor to. Since the amplitude of the saturation time is i? ~ 1 for correlated 
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models without additional crossover (BD, RDSR, RSOS, CRSOS, among others), we expect 
that 

B r^to (18) 

when the crossover from random growth is present. 

Now substituting the amplitude B from Eq. (|THj) in the Family- Vicsek relation (j21) and 
considering that f{x) ~ in the growth regime {to <^ t <^ t^), we obtain W ~ i^^^- 
Comparison with Eqs. (jH) and (fTTj) immediately leads to 

A ~ to'/'. (19) 

However, this relation may be obtained from different but consistent arguments, as follows. 
During time to the neighboring columns randomly grow, thus the local roughness Wi is of the 
order of the RD roughness to^^^. Wi represents height fluctuations within narrow windows, 
whose sizes are of the order of one lattice unit, in a large lattice and at times long enough 
for significant correlations to appear inside those windows. On the other hand, when the 
global roughness W attains saturation, the whole system is also highly correlated. Thus, if 
the total system size is L ~ 1 (i.e. system size is small but not equal to 1), this system is a 
narrow window, and we expect Wgat ~ Wi. Since Wsat ~ A in this case (Eq. El with L = 1), 
we obtain Eq. ^T^i . Certainly this argument does not apply to systems with anomalous 
scaling, where local and global roughness have different scaling properties, but that is not 
the case of the lattice models analyzed here or in related works. 

Now we consider the particular properties of the lattice models, focusing on their small 
length-scale features. Our arguments are similar to those of Braustein and Lam for 
the RD-BD and the RD-RDSR models, but here we will emphasize the independence of 
the results on the universality class of the process. In all cases, correlated growth attempts 
have probability p, while uncorrelated growth takes place with probability 1 — p. Thus, in a 
given column, the mean time interval between two successive depositions that buildup the 
correlations along the substrate is r = 1/p. 

Since BD involves lateral aggregation, a single deposition attempt following this model 
rules introduces height correlations between that column and the neighboring ones. This is 
illustrated by the deposition in column 3 of Fig. 2a: the large height difference from column 
4 is immediately suppressed by lateral aggregation. Consequently, in the RD-BD model, 
correlations between neighboring columns are buildup within an average time interval r. 
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Consequently, we expect to ~ t in any spatial dimension. Using Eqs. @, (fTHj) and 
(fTU|) . we obtain the exponents 5 = 1/2 and ?/ = 1 for the RD-BD model, while Eqs. (|7j) and 
(fT7|) provide exponents 7 dependent on /? and, consequently, dependent on the substrate 
dimension. 

Now we consider the RD-RSOS model as a typical example of solid-on-solid model. In this 
case, a single RSOS attempt is not sufficient to balance the random growth of neighboring 
columns. For instance, consider the aggregation rejection in column 3 of Fig. 2b: within 
the time interval r, a single particle is expected to be deposited at columns 2 and 4, but 
these events do not suppress the large height difference from column 3. Thus, within a 
time interval of order r, significant local correlations are not generated. Instead, in order 
to the rejection mechanism to balance the random growth locally, it is necessary that the 
number of rejected attempts at a given column be of the same order of the height difference 
between the neighbors. While the number of RSOS attempts at a given column during time 
t increases as tp, the local height difference in random growth increases as t^^"^. Matching 
these values we obtain the crossover time to ~ p^"^ for the RD-RSOS model. This result 
gives 6=1 and y = 2 in all dimensions. The exponent 7 depends on the substrate dimension 
due to its dependence on the scaling exponent /5 (Eqs. 171 and IT7|) . 

The same arguments apply to the crossover from random to EW scaling observed in the 
RD-RDSR model 0, Illl, since a single diffusion event to a lower height does not suppress 
a large height difference immediately. In the RD-CRSOS model, there is a combination of 
rejection of aggregation at high columns and diffusion to plateaus or valleys, but this com- 
bination also does not introduce significant height correlations immediately. Consequently, 
to ~ for these models, which lead to 5 = 1 and y = 2. 

The features of the generalized RSOS model are explained by extending the arguments 
ofRef. ^|. Chien et al [19| argued that the KPZ growth takes place when the randomly 
growing roughness (Eq. is of order S, which gives to ~ "S*^. This is the same form of the 
previous solid-on-solid models, with 5* interpreted as a reciprocal probability of rejecting the 
aggregation. Thus, 6=1 and ?/ = 2 for the generalized RSOS model. 

All the above results are in full agreement with simulation data shown in Sec. II. Although 
our simulation results (Sec. II) were limited to one-dimensional systems, Horowitz and 

n 

Albano [11| presented simulation results for the RD-BD and the RD-RDSR models in ci = 1, 
d = 2 and d = 3, an in all cases they agree with our predictions. In fact, no reference to 



11 



a particular system dimension was done in the above arguments. Instead, they were only 
based on random growth properties (Eq. P), the Family- Vicsek relation (Eq. ^ and the 
hypothesis of the existence of a crossover from random to correlated growth. 

Our analysis clearly separated solid-on-solid models with single particle aggregation at- 
tempts and models with lateral aggregation. Other limited mobility growth models may be 
classified in one of these two groups by inspection of the microscopic aggregation rules of 
the correlated process. Since the exponents 5 = 1/2 and y = 1 are found in models with 
some type of lateral growth and, consequently, excess velocity, they are typical of models 
in the KPZ class. However, models in 5 = 1 and y = 2 may be found in any universality 
class of interface growth, including KPZ. Moreover, although most systems previously stud- 
ied are two-component models, the application to the generalized RSOS model shows that 
the crossover exponents and the above scaling approach may be extended to other systems 
presenting crossover from random to correlated growth. 

Now we consider the particular case of KPZ systems in = 1, where relations between 

n 

scaling amplitudes and the coefficients of the growth equation are known : 

A ~ (20) 
5~z/i/2|A|-\ (21) 

and 

C~z/-'/=^|A|'/'. (22) 

In systems with a crossover from random to correlated growth, A and u may be arbitrarily 
small (but non-zero), while the noise amplitude D is finite, thus the dependence on D is 
omitted in Eqs. (HOI), (EH) and (HH). From Eqs. ^ and we obtain 

to ~ (23) 

and using Eqs. (fTHj) and (PT|) we obtain 

|A| ~ (24) 

in the crossover region of small u and small |A|. 

These relations may also be obtained from the condition that the random-KPZ crossover 



takes place at the same time of the EW-KPZ crossover, tc ~ (^^j ~ z^^A ^ 3, 



36 



37 
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In other words, after leaving the random growth regime, those systems do not show an 
additional crossover from EW to KPZ because the linear and nonlinear effects simultaneously 
turn up. 

Eqs. (I23j) and (j24j) provide relations between |A| and the model parameters p or S. 
u p and |A| ~ p^^"^ are obtained as particular relations for the RD-BD growth, but cannot 
be viewed as universal relations for the random to KPZ crossover m d = 1. It contrasts 
to what could be naively believed from some works on competitive models with random 
to correlated crossover Instead, a large variety of KPZ processes, such as the 

RD-RSOS and generalized RSOS models, follow the relations z/ ~ and |A| ~ p'^ in the 
crossover region m d = 1. At this point, it is important to recall that previous numerical 
work on the RD-BD model did not calculate v and A directly from the growth velocities and 

nn 

interface shapes [10|, [111 . Instead, their relations with the probability p were obtained from 
the scaling relations ()20|) . ()2H). and ()22j). similarly to what was done in other competitive 
models jl^ and in the present work. 



IV. CONCLUSION 



We studied limited mobility growth models with crossover from random to correlated 
growth. Universal relations between the crossover time to and the amplitudes of saturation 
roughness and saturation time were obtained from random deposition properties and Family- 
Vicsek scaling, for any spatial dimension. The lattice models with that crossover can be 
separated in two groups: the first one with lateral aggregation, such as ballistic deposition, 
in which correlations spread faster, and the second one of the solid-on-solid models with 
single particle aggregation attempts, which require much longer times for the correlated 
aggregation to balance the random growth. While to ~ P in the first group, where p is the 
small probability of the correlated mechanism to work, in the second class we showed that 
to ~ p^. These relations are independent of the universality class of the process: although the 
first group is expected to include only KPZ processes, due to the presence of lateral growth, 
several KPZ models are also found in the second group. All these features are confirmed 
by simulations of lattice models in three different universality classes and different substrate 
dimensions, partly obtained from other authors work. In = 1, relations between the 
probability p and the coefficients of the KPZ equation can be obtained for both groups of 
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models. 

Many natural and artificial processes involve the competition of different growth dynam- 
ics, and their analysis may eventually be improved by extensions of our scaling arguments. 
In the case of KPZ scaling, although the connection between scaling amplitudes and the 
coefficients of the continuous equation is not trivial in d > 1, the numerical calculation of u 
and A may help the search for those relations, at least in the small z/, small A limit. 
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FIG. 2: Aggregation rules of lattice models in d = 1, in which full squares represent aggregated 
particles, open squares represent incident particles, with the column of incidence indicated by ar- 
rows, small bullets indicate the aggregation position of the incident particles and crosses indicate 
rejected attempts of deposition. The rules of ballistic deposition are illustrated in (a). The condi- 
tions for RSOS aggregation (probability p) in the competitive model RD-RSOS are illustrated in 
(b). The generalized RSOS model with S = 3 is illustrated in (c). Column labels are indicated 
below the substrate lines. 
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FIG. 3: Time evolution of the roughness of the RD-RSOS model in lattices with L = 256, for two 
different probabilities p. The inset shows a linear plot of the same quantities. 
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FIG. 4: (a) Extrapolation of the ratio Wsat/L°' to L ^ oo, with a = 1/2, for the RD-RSOS model 
in d = 1 with p = 0.25 and 64 < L < 512. (b) Effective exponents 6p versus squared probability 
p"^ for the RD-RSOS model, (c) Extrapolation of the ratio W/t^ to t oo, with /3 = 1/3, for the 
RD-RSOS model in d = 1 with p = 0.2 and L = 8192. (d) Effective exponents 7p versus squared 
probability for the RD-RSOS model, with a least squares fit of the data. 
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FIG. 5: (a), (b): effective exponents 6s and 75 versus 1/S for the generalized RSOS model; (c), 
(d): effective exponents 5p and 7p versus squared probability for the RD-CRSOS model. The 
dashed lines in (a), (b) and (d) are least squares fits of the data. 
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